Model Independent Extraction of S-Matrix Poles from Experimental Data 
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By separating data points close to a resonance into intervals, and fitting all possible intervals 
to a simple pole with constant coherently added background, we obtained a substantial number of 
convergent fits. After a carefully chosen set of statistical constraints was imposed, we calculated the 
average of a resonance pole position from the statistically acceptable results. We used this method to 
find pole positions of Z and N(1440) resonances, and to show that the strong discrepancy between the 
old and new measurements of the T(11020) mass stems from specious comparison of the T(11020) 
pole with its Breit-Wigner mass. 
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Breit-Wigner (BW) parameters are often used for the 
description of unstable particles (see e.g. Review of Par- 
ticle Physics [l|), although shortcomings of such choice 
have been pointed out on numerous occasions. For ex- 
ample, Sirlin showed that the BW parameters of the Z 
boson were gauge dependent [2J]. To resolve this issue 
he redefined BW parameters, but also suggested usage of 
the S-matrix poles as an alternative, since poles are fun- 
damental properties of the S-matrix and therefore gauge 
independent by definition. In a somewhat different study, 
Holder advocated using S-matrix poles for characteriza- 
tion of nucleon resonances Q in order to reduce confusion 
that arises when different definitions of BW parameters 
are used [4|. However, loosely defined [5( BW parame- 
ters of mesons and baryons are still being extracted from 
experimental analyses, compared among themselves [l|, 
and used as input to QCD-inspired quark models [6| and 
as experiment-to-theory matching points for lattice QCD 
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The main motivation for this research are strong dis- 
crepancies between the old and newly obtained param- 
eters of some well known resonances. In particular, 
BaBar collaboration recently reported that the mass 
of T (11020) is 10996±2 MeV @], significantly different 
from the old value of 11019±2 MeV Q. Furthermore, 
the width turned out to be less than a half of the ex- 
pected, namely 37±3 MeV instead of 79±16 MeV. 

In this paper, we developed a reliable method for 
model-independent extraction of S-matrix pole positions 
directly from the data, and connected them to the Breit- 
Wigner parameters in order to understand the observed 
discrepancies. Wc showed that the both results, BaBar 
and PDG, are consistent: the old ones should be inter- 
preted as BW mass and width, while the new ones are 
pole parameters. 

The first step in devising a method for extraction of the 
pole parameters from the experimental data is to set up 



an appropriate parameterization. The parameterization 
presented here is based on the assumption that close to 
a resonance, the T matrix will be well described with 
a simple pole and a constant background. The similar 
assumption was used in Holder's speed plot technique [3[. 
The speed plot is a method used for the pole parameter 
extraction from the known scattering amplitudes. It is 
based on calculating the first order energy derivative of 
the scattering amplitude, with the key assumption that 
the first derivative of the background is negligible. 

The T matrix with a single pole and constant back- 
ground term is given by 
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where W is center-of-mass energy, r& and bb are complex, 
while M p and T p are real numbers. Total cross section is 
then given by a » \T\ 2 /q 2 , where q is the initial center-of- 
mass momentum. Equation ([T]), as well as other similar 
forms (see e.g. [1|), are standardly called Breit-Wigner 
paramcterizations, which can be somewhat misleading 
since M p and T p are generally not Breit-Wigner, but pole 
parameters (hence the index p). The square of the T 
matrix defined in Eq. ([1]) is given by 
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where, for convenience, we simplified the numerator by 
combining the old parameters into three new real-valued 
ones: X^, M z , and T z . Pole parameters M p and T p are 
retained in the denominator. 

With such a simple parameterization, it is crucial to 
use only data points close to the resonance peak. To 
avoid picking and choosing the appropriate data points 
by ourselves, we analyzed the data from a wider range 
around the resonance peak, and fitted localy the param- 
eterization ^ to each set of seven successive data points 
(seven data points is minimum for our five-parameter 
fit). Then we increased the number of data points in 
the sets to eight and fitted again. We continued increas- 
ing the number of data points in sets until we fitted the 



whole chosen range. Such procedure allowed different 
background term for each fit, which is much closer to re- 
ality than assuming a single constant background term 
for the whole chosen data set (see e.g. discussion on the 
problems with speed plot in Ref. [9[). In the end, we 
imposed a series of statistical constraints to all fits to 
distinguish the good ones. 

In order to pinpoint the statistical strategy to be used, 
we did a substantial number of simulations with the data 
sets that had known poles and zeros. It turned out that 
the most successful strategy was to make an ordered list 
of all fit results, from best to worst, and then to drop 
the worst three quarters using the followin g go odness-of- 
fit measures: Akaike information criterion [11], Schwartz 
(Bayesian information) criterion [l2|, and P-values of 
the extracted fit parameters (in particular, M p and T p ). 
Eventually, we kept the intersection of the fits that sat- 
isfied all criteria. Results closest to the original poles 
were produced by averaging the obtained pole positions 
of all good fits. The standard deviation turned out to be 
a good estimate for errors of obtained parameters. All 
other approaches we tested, such as keeping only a hand- 
ful of the best fits, or keeping just those whose values 
of reduced x 2 were close to one, failed to accurately re- 
produce the originalpole parameters. The whole analysis 
was done in Wolfram Mathematica 7 using Nonlinear- 
ModelFit routine [H]. 

Having defined the fitting strategy, we tested the 
method by applying it to the case of the Z boson. The 
data set we used is from the PDG compilation [1], and 
shown in Fig. [TJ Extracted pole masses are shown in the 
same figure: filled histogram bins show pole masses from 
the good fits, while the empty histogram bins are stacked 
to the solid ones to show masses obtained in the discarded 
fits. Height of the pole- mass histogram in Fig. [TJ is scaled 
for convenience. 

Extracted S-matrix pole mass and width of Z boson 
are given in Table HI The pole masses are in excelent 
agreement, while the pole widths are reasonably close. 
It is important to stress that the difference between the 
pole and BW mass of the Z boson is fundamental and sta- 
tistically significant. Distribution of discarded and good 
results is shown in the lower part of Fig. [TJ 



TABLE I: Pole parameters of Z obtained in this work. PDG 
values of pole and BW parameters are given for comparison. 



Pole Pole PDG [1] BW PDG [1] 



M/MeV 
T/MeV 



91159 ± 8 
2484 ± 10 



91162 ± 2 
2494 ± 2 



91188 ± 2 
2495 ± 2 



Next, we turn to the data from BaBar collaboration 
8] to determine whether the PDG averages for T(11020), 
or the newly reported resonance parameters obtained in 
Ref. [al are correct. 
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FIG. 1: [Upper figure] PDG compilation of Z data [l[ and 
histogram of obtained pole masses. Line is the fit result with 
the lowest reduced \ 2 (just for illustration). Dark (red online) 
colored histogram bins are filled with statistically preferred 
results. [Lower figure] Pole masses vs. pole widths. Dark 
(red online) circles show statistically preferred results we use 
for averages. 



Our local fit of the T(11020) pole parameters is shown 
in Fig. [21 As in the case of Z boson, the full and empty 
histograms show how many of the extracted pole mass 
fits were accepted or discarded in the analysis. 

Average values obtained for the resonance mass and 
the width are given in Table [Hi together with the same 
parameters obtained in the BaBar analysis, and those 
quoted by PDG. Extracted pole parameters of T(11020) 
are practically the same as those reported in |8j, even 
though our parameterization is much simpler (single pole 
plus constant background vs. two constant background, 
and two pole terms). The ori gina l results cited in PDG 
(from CUSB []]|, and CLEO [11 ) were obtained by fit- 
ting Gaussians to the resonance peaks in the data, and 
peak positions are usually closer to the BW mass. 

To investigate this case further, we analyze another 
resonance with strong difference between the pole and 
BW mass, the Roper resonance N(1440). We extracted 
Roper resonance pole parameters from the ttN elastic P\\ 
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FIG. 2: T(11020) resonance pole obtained by our method. 
Grey (green online) rectangles represents PDG range for 
T( 11020) mass (both figures) and width (only lower figure). 



TABLE II: Parameters of T(11020) meson. Pole parameters 
are results of this work. 



T(11020) 



Pole BaBar [1, 8] 



M/MeV 
T/MeV 



10999 ± 1 
38 ± 1 



10996 ± 2 
37 ± 3 



PDG [1] 

11019 ± 8 
79 ± 16 



TABLE III: N(1440) resonance parameters. 



iV(1440) 



Pole Pole PDG [1] BW PDG [1] 



M/MeV 
T/MeV 



1370 ± 6 
197 ± 6 



1365 ± 15 
190 ± 30 



1440 ± Ig 
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FIG. 3: N(1440) resonance pole obtained by our method. 
Grey (green online) rectangles represents PDG range for 
N(1440) Breit-Wigner mass (both figures) and width (only 
lower figure). 



partial wave obtained in the GWU analysis [10( . Accord- 
ing to PDG, this wave has a very rich structure: there 
is a four-star Roper resonance, a three-star N(1710) res- 
onance, and a one-star N(2100) resonance. However, the 
GWU analysis reports only one resonance in this par- 
tial wave, the Roper N(1440) resonance. In a prelimi- 
nary analysis, we could see some indication for all reso- 
nances mentioned by PDG but, for this study, we focus 
on N(1440) because of its unusually strong shift between 
the pole and BW mass (roughly 75 MeV). 

Our results for N(1440) are given in Table Mil where 
we see that the pole parameters are in an excellent agree- 
ment with the PDG estimates. Unlike the pole mass, BW 
masses are situated closer to the positions of the peak (see 
Figs. H and U). 

The field-theory reason for the resonance pole shift is 



the energy dependence of the imaginary part of resonance 
self energy, which is commonly modeled by the energy 
dependent width [18|, |l9( . Equation ((T|) would be exact if 
the self energy was constant. In more realistic cases, the 
T-matrix denominator D(W) is given by 



D{W) =M b -W -iT(W)/2, 



(3) 



where we introduce the BW mass and Mb, which is gen- 
erally not the real part of the pole position. Keeping only 
the first two terms in Taylor expansion of r(l4 7 ) about 
W = Mb (as done in Ref. [19|) the width becomes 



T(W)/2 = T b /2 + t&nd (W -M b ), 



(4) 



where the T b is a shorthand for T(Mb), and tan# is the 
slope of r/2 at W = M b . The pole position M p - iT p /2 



is obtained by solving D(W) = 0, which yields 

AI P = M b + sin 9 cos 9 T b /2, (5) 

r p = cos 2 6»rV (6) 

Relations ([3]) and © , originally introduced in Ref. [19f , 
may be used to cross check estimates for pole and BW pa- 
rameters. Table iPVl shows angles 9 for all resonances an- 
alyzed in this paper, calculated from PDG estimates for 
pole and BW widths using Eq. ([6]). Since M p is smaller 
than Mb, we chose negative 9 solution (cf. Eq. ([5]). It 
turns out that the BaBar value of Y(11020) mass is 
accurately reproduced. 



We are today witnessing the dawn of ab-initio calcu- 
lations in low-energy QCD. In order to compare theo- 
retical predictions with experimentally determined res- 
onance states, we need first to establish proper point of 
comparison. The case of T(11020) is a vivid example how 
particularly careful we must be when choosing this point. 
Therefore, we would like to express our concern about 
other potentially problematic comparisons between the 
pole and BW parameters in the literature, in particular 
in the Review of Particle Physics, and recommend draw- 
ing a clear distinction between the two in future publica- 
tions. 



TABLE IV: The connection between S-matrix pole and Breit- 
Wigner parameters using only the PDG values. 





e/° 


M, 


,/MeV PDG|JJ 


M p /MeV Eq. © 


A(1232) 


-23.0 




1210 ± 1 


1210 


A/"(1440) 


-37.3 




1365 ± 15 


1368 


T(11020) 


-46.8 




10996" ± 2 


10999 


Z 


-1.26 




91162 ± 2 


91161 



"BaBar value. 
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Does this 9 carry any physical meaning? For res- 
onances with one dominant decay channel, such as 
the A(1232), we can impose a unitarity condition 



(ImT 



\T\ 2 ) to Eq. flU and learn that 



■no 



and b p = e sin 9. It is the same 9 and represents a half 
of the complex residue phase. Indeed, from Ref. [l[ we 
read that A(1232) has (— 47 ± 1)° for pole residue phase, 
quite consistent with -46°, a double value of the 9 from 
Table ITVl However, this simple relation is lost when im- 
portant inelastic channels are open, e.g. in the N(1440) 
case, where 2 9 « -75°, which is significantly larger than 
its residue phase -100° [l|. The difference comes from 
different r(W) in the denominator and numerator of T 
matrix: total decay width is in the denominator, while 
the partial widths are in the numerator. Since the energy 
dependence of the two is in general different, their slopes 
(i.e. tan 9) will be different as well. 

Since our pole extraction method confirmed BaBar 
result, the successful cross check is the last piece of the 
puzzle. PDG estimates of T(11020) are consistent with 
BW parameters. 

In conclusion, we have developed a model-independent 
method for extraction of resonance pole parameters from 
total cross sections and partial waves. Very good esti- 
mates for Z boson, T(11020), and N(1440) pole positions 
were obtained. Furthermore, we showed that the strong 
discrepancy between PDG estimates and BaBar result 
for T(11020) comes from specious comparison of the pole 
and BW mass. 
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